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Abstract 

In this paper we examine the supersymmetric Lee -Yang model in the 
presence of boundaries. We determine the reflection factors for the Neveu- 
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1 Introduction 



Boundary conformal field theories attracted much interest recently, due to their 
applications in describing D-branes in string theory [U [2] and their relevance in 
condensed matter physics, e. g. in the Kondo problem [3]. 

Several papers appeared in the literature that deal with the consistent boundary 
conditions and reflection factors pQ EI El d El [9] , the boundary perturbations and 
the corresponding renormalisation flows pm [EU E21 CE] of two-dimensional super- 
symmetric field theories. 

One of the simplest supersymmetric theory is the supersymmetric Lee -Yang 
(SLY) model and in this paper we study this model in the presence of boundaries. 
In a paper [12J Ahn and Nepomechie proposed reflection factors for the supercon- 
formal boundary conditions and determined a boundary flow using Boundary Ther- 
modynamic Bethe Ansatz (BTBA). Here we argue that both their reflection factor 
and one of their TBA formulae are mistaken. We propose different reflection factors 
which we check using Boundary Truncated Conformal Space Approach (BTCSA). 
We also predict a different boundary flow using BTBA and we compare the fixed 
points and the change of the boundary entropy again with our BTCSA results. 

The paper is organized as follows. In section 2 we summarize some basic facts 
about the SLY model. In section 3 we determine the reflection factors by applying 
a "folding" trick to the supersymmetric sine-Gordon reflection factors. In section 4 
some details of the BTCSA are explained briefly and after calculating the massless 
reflection factors the predicted energy levels and phase shifts are compared to the 
BTCSA results. In section 5 we write down the BTBA equations and calculate the 
variation of the ^-function along the boundary flow. In section 6 this is compared 
to the BTCSA result for the fixed points and the g-function. Finally, in section 
7 we extend our study of the flows to the models SM(2,kn + 4) which are the 
generalisations of the SLY model. 

2 The supersymmetric Lee -Yang model 

The supersymmetric Lee -Yang model or the superconformal minimal model 
SM(2,8) is a non-unitary conformal field theory with central charge c = — The 
super Kac-table with the highest weights of the primary fields is 
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The algebra of superconformal transformations in the plane are generated by two 
fields, T(z) and G(z) whose mode expansions read 




(2.1a) 
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The N = 1 supersymmetric extension of the Virasoro algebra is defined by the 



following (anti) commutation relations: 

[L n , L m ] = (n - m)L n+m + ^n{n 2 - l)5 n _ m , (2.2a) 

{G r , G s } = 2L r+s + ^{r 2 - 1/4) <5 r ,_ s , (2.2b) 

[L n , G r ] = (~-r) G n+r . (2.2c) 



Here n, m denote integers, while depending on the boundary conditions for G(z) 
the indices r, s can take half-integer (Neveu-Schwarz sector) or integer (Ramond 
sector) values. 

The highest weight representations of the algebra are defined by highest weight 
states | A) satisfying 

L n |A)=0 n>0, (2.3a) 

L |A) = A, (2.3b) 
G r |A) = r>0. (2.3c) 

The Verma module is generated by operators having a nonpositive index. The 
vectors 

L ni ...L nk G n ...G n \A) , ni<---<n fc <0, n < • • • < n < (2.4) 

constitute a basis for the Verma module. Note that the inequalities for the r^-s are 
strict since Gf = L 2r — ^ <5 r ,o by equation (I2.2bl) . 

In this paper we consider the so-called spin model which can be obtained by the 
projection onto the states with even fermion parity of the Hilbert space. 

One can consider this superminimal model on a strip of width R with nontrivial 
boundary conditions at the edges of the strip. In the Ramond sector the superconfor- 
mal boundary conditions are in one-to-one correspondence with the highest weight 
representations so they can be indexed by the same label set (the Kac indices), 
whereas in the Neveu-Schwarz sector there are two boundary conditions for each 
entry of the Kac table, the NS and the NS boundary conditions [4]. 

This model can also be defined as the minimal model M(3,8) with Kac table 
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Suitable combinations of the Virasoro representations constitute irreducible repre- 
sentations of the superconformal algebra: 

_(-!) ffi(i) (2.5b) 

V 4/SVir V 4/Vir V4/Vir V ' 

^\ <*(™\ . (2.5c) 



32/SVir V 32/Vir \32/Vir 

<—(--) • (2.5d) 

32/SVir V 32/Vir v 7 

Since the superprimary fieldfare also Virasoro primary fields and the g-factors are 
the same for the NS and NS boundaries, the ^-factors of the various conformal 
boundary conditions can be calculated in the Virasoro picture. The matrix elements 
of the S modular transformation are 

S( r ,s),(r>,s>) = 2^-(-l) rs ' +r ' s+1 sin [rr'%^\ sin ( 2 - 6 ) 

and according to [T4] the g-factors are given by 

V Pn.ol 

where denotes the conformal vacuum and fl denotes the state of lowest conformal 
weight. In our case = (1, 1) and Q = (1, 3), so 

(-l)* sinfs 

If we choose the (3 = (1, 1) boundary condition at one of the edges and a = (1, s) 
at the other one, the partition function on a cylinder of circumference L and length 
R will be 

Z{R,L) =tr NS 1 + ( ~ 1)F e- LH (^) + tr R 1 + V lh (m),(m) = 



l±f-l 



tr (ljS) ^_ e -*sn*-<v"; , (2.9) 

where the sign is "+" if both boundaries are of NS or NS type and "— " if they are 
different (see [HHH] for details). Thus the Hilbert space consists of the bosonic levels 
of the corresponding superconformal Verma module if the boundary conditions are 
both of NS or NS type and of the fermionic levels of this module if they are different. 
For a bosonic primary field the bosonic levels are the integer levels and for a fermionic 
primary field they are the half-integer levels. 
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We consider the boundary perturbation of this model defined by the action: 

oo 

A = A SM{2 , 8 ) + h J dt(G_ 1/2 li3 )(O,t). (2.10) 

— oo 

where (G_i/20) is the field appearing in the operator product expansion of G(z) and 
0: 

m w \ (^-1/20)0) ( . 

G(z)(p(w) = '- h • • ■ (2-11) 

This operator can live on the (1, s) boundary since the fusion coefficient iVl^w^j 
is non-zero. This boundary perturbation is relevant, since the conformal dimension 
of the perturbing operator is 1/4 < 1. Since there is no perturbation in the bulk, 
the theory remains massless along the flow. The perturbation preserves supersym- 
metry so the boundary renormalisation group flows will take place in the space 
of possible supersymmetric boundary conditions, in which the fixed points are the 
superconformal ones. This perturbation is also integrable [T5] . 



3 The reflection factor 

We obtain the S-matrix and the reflection factor of the supersymmetric Lee -Yang 
model by considering it as a reduction of the supersymmetric sine-Gordon (SSG) 
model. The original idea is of Smirnov [16j and it was first applied in supersymmetric 
context for bulk S-matrices in [PT] . The SSG S-matrix was also derived in [T7], the 
formulae are collected in appendix lAl 

If A > 2 in the sine-Gordon model the pole in the S-matrix describing the 
scattering of the first breather on itself corresponds to the second breather as a 
bound state. For A < 2 only the first breather is in the spectrum and the pole 
is explained by the Coleman-Thun mechanism |18| which requires the presence of 
solitons in the theory. However, at the particular value of the coupling constant 
A = | the masses and the S-matrices for the first and the would-be second breathers 
become equal: mi = m-i and Sl' SG = S$$ G = S^ G , so the pole can be explained by 
the self fusion of the only breather and the solitons can be projected out consistently 
from the theory [T7]. Thus the SSG model is reduced to the supersymmetric scaling 
Lee -Yang model containing only one superdoublet. The S-matrix simplifies to 

W*) = S ly (0)Ssvsy(0) (3.1) 

where 



sinh(#) +isin( 2Tf 
sinh(6 l ) — i sin( 



SlyW = -TTT^ • (3-2) 
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and 



smh(6) 



SUSY 



(6) = Z{6) x 



sin(J_) 



Z{9) 




sin(f) 
cosh(|) 







-1 + 2% - 



sin(f) 



MO) ) 



"2^ 



sinh(|) +zsin(f 



dt sinh(g) sinh( 
t cosh 2 (4)cosh(t) 



smh( — J 

7T 



(3.3) 



(3.4) 



This is the S-matrix used also by Ahn and Nepomechie in [12]. For later convenience 
we introduce the notation 

Z(9) = S hY (8)Z{6) . (3.5) 

It is easy to check that this S-matrix satisfies the appropriate unitarity and 
crossing relations and that the Yang-Baxter equations are also satisfied. 



3.1 The "folding" 

We have seen that at the particular coupling A = § the bulk SSG model reduces 
to the SLY model. For this reduction to be consistent also in the boundary case it is 
necessary that at this point the reflection factors for the 1st and the would-be 2nd 
breather be equal. This requirement gives functional relations between the two SSG 
boundary parameters, rj and $ (first considered in [H]). 

The reflection factors of the boundary SSG model can be found in appendix [Bl 
In the BSSG + case the supersymmetric factor -Rsusy(#) (see IB.8I) does not contain 



the boundary parameters and the only way it depends on the species of the breather 
is through the parameter p which at this special coupling is pu = vr — k^f. From 
pi = —p2 = | it is obvious that A±{6) are automatically the same for k = 1, 2. 

Now we turn to the scalar part, i. e. to the SG reflection factors. It turns out 
that RW(6) = RW(9) if 

2k + 1 

v = i$ + — - — 7T, keZ. (3.6) 

If k = 1, i. e. rj = {■& + ^ then = = R where 

where we introduced the notation 



2/ V2 



sinf tt: + 
sin(|-x|; 
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If k = 1, i. e. 7] = i-& + | then = #( 2 ) = i? 7 with 
l + 6\ /l-&\ -1 /5-&\ f5 + b 



R T (6) = R (6) 



2 



M0)[s LY {e + t b -^n)s LY (e-z b -^^ , (3.8) 



where we introduced the parameter b = — 2. Rq{6) and Ri{6) are the same 
reflection factors as those given by Dorey et al. (for example in [T9]). 
Finally, if k — — 1, i, e. r] — i$ — f then .R^ = i?( 2 ) = R n where 



i? 7J (0) = i2o(0) 



6-l\ /6 + 3\ 1 fb-3\ /5 + b\ 1 



^)^Y(e + i^7r)5 LY (e-i^l7r)^ . (3.9) 



The other values of k do not lead to new reflection factors. With the redefinition 
b ^ b + 4 the factor Rn{6) becomes 



R n {9) = Ro{6) 



b + 3\ f5-b\ fb + l\ (3-b 



Ro(e)(s LY (e + z b ^7ry LY (e-i b -^^ , (3.10) 



which can be identified as the reflection factor for the excited boundary [19]. Thus Ro 
and Rj are the two ground state reflection factors for the two boundary conditions 
of the Lee -Yang model. They satisfy the boundary unitarity and boundary crossing 
equations: 

R(0)R(-G) = 1, (3.11) 

r (t~°) = s M20)r(y + 6 )- ( 3 - 12 ) 

The SLY reflection factors also satisfy the boundary unitarity equation provided 
A±{6)A±{—6) = 1 and this is indeed true: 

A ± (9)A ± (-9) - 2 — — 1 . (3.13) 

The boundary crossing equations are 

2"-^(f -e) = s**{20)K(t + e ) + 5 » 2 *K(y + d ) > ( 3 - 14a ) 
2"-^(f -e)= s$}(20)H*(^ + e)+ s**(20)b*(^. + e) . ( 3 .i4b) 
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Since the Lee -Yang parts satisfy the scalar version of the boundary crossing condi- 
tion we only have to check that the supersymmetric factors Ssusy and -Rsusy satisfy 
(I3.14p . Using the integral representations for them this can be proved analytically. 
We also checked that the matrix parts satisfy the boundary Yang-Baxter equations 
as well. 

It seems plausible that the two supersymmetric reflection factors, R (9)Rs\jsy(9) 
and Ri (9) Rsusy (9) are the reflection factors for the two Neveu-Schwarz boundary 
conditions of the supersymmetric scaling Lee -Yang model, since for the Ramond 
type boundary conditions we do not expect diagonal reflection factors (the fermion 
parity can change). It is also natural to think that i? -RsusY is for the (1,1) boundary 
condition, because it does not depend on any boundary parameter. 

We emphasize that the reflection factors proposed here are different from the 
non-diagonal one proposed by Ahn and Nepomechie in |12j . 

3.2 Relation between the boundary parameter b and the 
boundary coupling h 

Following reference [20j it is possible to conjecture the relation between b and the 
coupling constant of the boundary perturbation, h (a similar argument was used also 
in [2T]). In that paper it was found that the transformation b — ► 4 — b is a symmetry 
of the model, since this change of b maps the reflection factors of the ground state 
and excited boundaries into each other. For b > 2 the pole at 9 = l ( b+ ^ n j n Rj 
which corresponds to an excited boundary leaves the physical strip but at the same 
time the pole at 9 = ^ 5 ~ b - )7r in Rjj enters. Thus for b > 2 R n describes the ground 
state boundary and Ri the excited boundary. A similar phenomenon is probably to 
occur in our case, but exactly the same thing can not work, since now the excited 
boundary state is a doublet while the boundary ground state is a singlet and Rn 
is a 4x4 matrix. However, at the pole Rn becomes a 1-dimensional projector. Since 
the reflection factors are manifestly periodic in b with a period of 12 there is a point 
(b = 8) beyond which the boundary energies swap back and the ground state will 
be a singlet again. This gives a symmetry b — > 10 — 6. Together with the trivial 
symmetry b — > —6 — b of the reflection factors this results in a periodicity of the 
physical spectrum in b with a period of 16. Similarly to the case in |20j , the value 
at half way between the fixed points of the two symmetries should correspond to 
h = 0. This leads to the conjecture for the relation between b and h: 



with h c = am 3 / 4 where a is some numeric constant and m is some mass scale of the 
model, e. g. the mass of the particle. 




(3.15) 
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3.3 The massless limit 

It is the supersymmetric part -Rsusy for which it is easier to calculate the massless 
limit since it does not contain boundary parameters. All we have to do is to take 
the limit 9 — ► oo. 

A±{9) can be written in the following form: 

A ± {6) = (3.16) 
y cosh(0) 



where 



with 



oo oo 



1(9) = I ^/(t)sinh(^) = -\J je^f(tn) (3.17) 



o 



Now 



so 



1 cosh(§) + cosh 2 (|) 
t[t) ~~l cosh 2 (|)cosh 2 (|) ' (3 - 18) 

Z7T i 

Aoo) = -y/(0) = -m (3.19) 
lim 2&A±(0) = V2 tl£l e x = e W) , ( 3 .20) 

which means that 

^susy(^) - (J fj ■ (3.21) 

Let us turn to the scalar parts. Ro(9) does not depend on the boundary param- 
eters and 

lim Ro(9) = 1 . (3.22) 

9— >oo 

Ri(9) does depend on the boundary parameter b. In the massless limit we want to 
keep h fixed, which means that the sine in (13.151) has to go to infinity in a proper 
way. This can be achieved only if b becomes complex. We should set b = — 3 — ib, 
then 

b'K . 3 H c jrw 3 

h = — /i c cosh( — )m 4 — ► (e em)*. (3.23) 

8 2 

This means that me^ should be kept fixed while m — > and 6 — > oo. We also 
want to keep the physical momentum, mcosh(0) fixed, which implies that me 9 is 
constant. We conclude that while 9, b — ► oo the combination — ~7r is fixed. So 

i?i(0) = R (9) (s LY (e + ^)^ly(^ - |tt)^ — Sg(0 - ^) = &9 " ^b) 

(3.24) 
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where 9 is the pseudo-rapidity and 9 B is related to 6.0 

We have found that the massless reflection factors for the two NS boundary 
conditions of the SLY model are 

R Q = (I fj and R T (9) = S h ^(9 - 9 B )R • (3.25) 
In the following we will omit the hat in the notation for the pseudo-rapidity. 



4 TCSA 

The method we use for examining boundary flows is the so-called truncated con- 
formal space approach, or TCSA. In this approach the infinite dimensional Hilbert 
space is truncated to a finite dimensional vector space by using only those states 
whose energy is not greater than a threshold value, E cut . This is equivalent to trun- 
cating the Hilbert space at a given level. The Hamiltonian is then diagonalised on 
this truncated space. The original idea was proposed in [22j, it was applied for the 
first time for boundary problems in |20j and in a supersymmetric theory in [23J. 
The current method is based on the techniques of [24], modified for superconformal 
minimal models for the first time in |13j . 

The Hamiltonian of the perturbed superconformal field theory takes the form 

H = H C ft + H peit , (4.1) 
where Hcft is the Hamiltonian of the unperturbed theory 

#cft = H a0 = ^(L - ±) (4.2) 

and .f/pert comes from the perturbation on the strip: 

fljf = MGU/ 2 0i, 3 )(O,O). (4.3) 

The location of the left boundary is at x = and we are free to choose t — for 
calculating the spectrum. By the exponential map this on the z-plane becomes 

tfpert = h (-) (G-l/201,3) (z = 1) . (4.4) 



Thus the complete Hamiltonian on the plane can be written as 

H = ^ L o-^ + h r~) (GU/2<M(1) ■ (4.5) 



lr The method used in [12], taking the joint limits 9 — > oo, m —> and demanding that the result 
be finite and unitary, leads to the same result. 
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Since a numerical calculation requires dimensionless quantities we have to introduce 
some mass scale, \i and measure the volume (R) and the energies in units of \i. In 
other words, we use the dimensionless quantities r = e = E//J,, n = /i//i 1//2 ~ Al - 3 
and h = H / ji: 



h 



71 

r 



c f r \ V 2- 3 
Lo ~24 + K U (G-i/ 2 0i, 3 )(l) 



(4.6) 



where k characterizes the relation between the boundary coupling and the chosen 
mass scale, thus it is not fixed to any preferred value yet. 

In order to calculate the eigenvalues of the matrix h we choose the basis vectors 
for the truncated Hilbert space of the form written in (12.40 with A = A ljS . After 
eliminating the null vectors and the corresponding submodules, we obtain a non- 
orthonormal basis, {e«}. Then the matrix elements of the TCSA Hamiltonian are 



h 



71 

r 



c / t \ 172-^1,3 



71 



where M is the inner product matrix 

Mij = (ei\ej) , 

and B contains the matrix elements of the perturbing operator: 



B, 



(4.7) 



(4.8) 



(4.9) 



These matrix elements can be calculated using contour integration techniques and 
the superconformal operator algebra, as worked out in detail in [1 3J . 



4.1 TCSA fits 

The Bethe-Yang equation for a massless particle moving between the boundaries 

e^RoRriO) = 1 , (4.10) 



is 



that is 



e^ R S^(6-6 B )R 2 (6) = l. 



Taking the logarithm we get 



where 



2pR + 5(9 - 9 B ) = 2n% 



8(9) = - log S£(9) 



(4.11) 

(4.12) 
(4.13) 



and n can be integer or half-integer for the two eigenvalues of Rq (1 and —1). 

For a massless particle p = E = ^e 9 where fi is some mass scale. A shift in the 
pseudo-rapidity 9 is equivalent with a change in \i. The shift in 9, 9 — 9 B can be 
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compensated by changing /i which for the TCSA means changing the numeric value 
of the dimensionless coupling k. 

We should note that the Bethe-Yang equations are not derived rigorously for 
massless particles. However, experience shows that it works and gives the correct 
energy levels even for rather small volumes. The massless limit of the massive TBA 
equations are the same as those derived from the massless Bethe-Yang equation and 
they work in a lot of cases, which is another reassuring fact. We will see that also 
in our case it gives consistent energy levels. 

Equation (14.121) can be solved for r (pR = ER = er) for different values of n and 
we can plot the inverse of the function r(e) and compare the lines with the TCSA 
data, £tcsa(t")- For this we have to distinguish the one-particle energy levels. For 
the multi-particle states we have a coupled system of Bethe-Yang equations that 
contain also the bulk S-matrix, which makes the behaviour of the energy eigenvalues 
different for finite volumes. This, in principle, allows us to select the one-particle 
levels. Furthermore, for large volume (r) the interaction becomes negligible and the 
asymptotic behaviour of the energies are 

E n ~ -n (4.14) 



for one-particle states and 

E ni ,n 2> ... ~ ~J2 Ui ( 415 ) 



for multi-particle states. Thus the IR behaviour is the same for one-particle and 
multi-particle states and for a given (half-)integer number n there is precisely one 
one-particle state and possibly several multi-particle states with asymptotic energy 
nix jr. If for a certain value of n there is only one level, one can be sure that it is 
a one-particle level. For n = \ there are no multi-particle states obviously and for 
greater values of n certain multi-particle states may be forbidden by exclusion rules. 
For example in the n — 1 case the two-particle state corresponding to 1 = | + \ is 
excluded, so for n — 1/2, 1 there is only one energy level which corresponds to a 
one-particle state. 

For each value of k first we must find the appropriate value of (9 B , for example by 
fitting the lowest energy level which is a one particle level. After this is done, all the 
other one particle energy eigenvalues should automatically fit the other lines. The 



result of such a fit can be seen in Figure 1(a) for the NS-NS case with k > 0, when 



only the integer levels of the module (1,3) are in the Hilbert space (see (12.91) and 
below). As can be seen, the TCSA spectrum possesses the expected features: the 
lowest energy levels are non-degenerate in the IR and the higher levels are arranged 
in groups. 

In every group there is exactly one level that fits the one-particle Bethe-Yang 
energy, with integer values n. It is interesting to observe that these levels are always 
the highest ones in their group, which is due to the fact that the interaction is 



attractive. The lines of Figure 1(a) not fitted by a solid line correspond to multi- 



particle states. The number of these levels is neither consistent with bosonic nor 
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with fermionic exclusion statistics. This suggests that these particles obey some 
generalised exclusion statistics, similarly to the ordinary Lee -Yang model [25]. 

One would expect that the other eigenvalue of the reflection factor corresponds to 
the reflection of the other memberof the superdoublet, that is the energy eigenvalues 
from the half-integer levels (NS-NS case) will fit the lines of half-integer values of 
n. Surprisingly, this is not the case and it is not clear why it does not happen and 
then what the reflection factor is for these fermionic levels. 

For k < (h < 0) the spectrum becomes complex, just like for the ordinary 
Lee -Yang model [20] . 

Another and more precise approach is obtaining the phase shift 5(6) from the 
various TCSA eigenvalues and comparing it with the exact function (14. 1 3B . From 
( 14.121) one gets 



<5TCSA(k)g 9 B ) 
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^TCsA(loge T csA - #b) = 2mr -2er. 



(4.16) 



It can be seen that our ignorance of the proper mass scale (or the value of 9b) is 
only a matter of a horizontal shift of the function along the #-axis while changing 



the quantum number n shifts the function vertically. In Figure 1 (b) the phase shifts 



calculated from the lowest one particle eigenvalues are plotted onto each other, which 
shows that the phase shifts extracted from the different one-particle levels are the 
same. This shows that the massless Bethe-Yang equation is meaningful and gives 
consistent results. The difference between the phase shifts for small rapidities is due 
to the truncation errors of the TCSA in the IR. 

In Tabled] the phase shift values calculated from the difference of the two lowest 
TCSA eigenvalues are compared with the theoretical values (14.131) . In Figures 



1 c 



and 1(d) the the TCSA phase shifts obtained from the lowest eigenvalues are plot- 
ted together with the proposed phase shift (14 .131) . As it can be seen there is good 
agreement. 





(a) Energy eigenvalues from the integer levels (b) <5tcsa(#) phase shifts calculated from dif- 

with the theoretical energies from the reflection ferent TCSA eigenvalues 

factor 
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#TCSA 


8(0) 


^tcsa(^) 


0.0025 


3.1358 


3.0401 


0.0570 


3.0099 


2.9276 


0.1154 


2.8760 


2.8076 


0.1782 


2.7337 


2.6796 


0.2462 


2.5824 


2.5426 


0.3204 


2.4217 


2.3961 


0.4022 


2.2513 


2.2391 


0.4935 


2.0707 


2.0707 


0.5968 


1.8795 


1.8898 


0.7158 


1.6772 


1.6953 


0.8562 


1.4630 


1.4860 


1.0275 


1.2361 


1.2606 


1.2461 


0.9951 


1.0179 


1.5461 


0.7378 


0.7567 


2.0163 


0.4612 


0.4765 


3.0696 


0.1609 


0.1765 



Table 1: The theoretical phase shift ([4T3l and the TCSA result 




0.5 1 1.5 2 2.5 3 3.5 4 0.5 1 1.5 2 2.5 3 3.5 4 



(c) TCSA phase shift from the difference of the (d) TCSA phase shift from the difference of the 
first two eigenvalues 3rd and 1st eigenvalues 

Figure 1: Reflection factor and TCSA spectrum 

5 TBA 

The Thermodynamical Bethe Ansatz is based on the quantisation of the energy 
levels using the S-matrices and the reflection factors of the model. The correspond- 
ing Bethe- Yang equation takes the following form in terms of the densities for N 
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particles in the periodic case ([26, 27J): 



2nP(9) = mcosh(9) + J d9' 



(5.] 

where p{9) denote the density of the occupied states and P{9) the density of available 
states. Similarly, po(9) denotes the density distribution of real solutions Xj of the 
equation 



in ^ 



n V tanh(^2^" ~ -qJ 
M tanh(^ + |)=- F ' (5 ' 2) 

(F is the fermion parity) for which €j = +1 and by po{6) the solutions for which 

9 = -I- 

The kernel functions are 

3. . / Z(») \ 



M9) = _I m l„ g ^_^J, ( 5 . 3a) 

where Z(#) is the scalar part of Ssly(#) defined in equation (13.51) . 

Let us turn to the boundary case. Since the reflection factors are diagonal they 
do not affect the calculation of the transfer matrix eigenvalues, so in the massless 
case 



p{9')$ z {9-9') 

+ Pomfae - e 1 ) + p {e'){- l -me - e') 



(5.4) 



where 



*(9) = ^lmlog(±R Mt (9 - 9 B )R light (9 - 9 B )) . (5.5) 
For Po(9) = po(9) + po(9) we have a similar equation, obtained from equation (15.21) : 

d T , /tanh(|(0-0'-if))' 



9' + ^)), 

' d9'p{9)${9-9'). (5.6) 



2rfoW = / deVW | Imlog (^M| 



Now using the fact that Pq(9) = po(9) + po(9) and equation (15.61) we can eliminate 
Po(9) from equation (15.41) getting 

m = + SI + ( p * ( ** ~ \* * $) ) {e) ~ {p ° * $) w ' (5 - 7) 
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where we introduced the notation 

/J/V 
— <f>{6 - e')W) (5.8) 

for the convolution. It turns out that the following identity holds ([12]): 

($,-!$*$)(#) = $(#). (5.9) 
The free energy of the system is given by 

f = J d6{p(6)ve d 

-T[P{6) log P(9) - p(6) logp(fl) - (P(9) - P {6)) \og{P{9) - p(0))] 

- T[P (8) \ogP (9) - po(9) logpo(^) - (Po(9) - p (9)) \og(P (6) - p (6))] } . 

(5.10) 

Now the densities p(9), po{9) can be varied to minimize the free energy. Using 

8P = 5p*<S>-8p- *$, (5.11) 
5P = -5p*§ (5.12) 



and introducing the quasi-particle energies 

p{9) e~ e W p (9) e- e °W 



P{9) 1 + e~< 9 ) ' P o (0) 1 + e- £ oW 

we arrive at the following TBA equations: 



(5.13) 



e(8) = ^e e -(<S>*(L-L ))(8), (5.14a) 
e o (0) = ($*L)(0), (5.14b) 

where L{9) = log(l + e~ e ^), L (9) = log(l + e~ e °^). 

Now using these equations we can write the extremum of the free energy: 

- = [ d9(pe e R + V(8))L(6) . (5.15) 



T 2vr _ 

Since the partition function on the cylinder behaves for large R as 

log Z aP = ~ » log(g a g p ) - RE™ = log(g a9p ) (5.16) 

we obtain 

oo 

log(^) = ^ J d8 (0) log(l + e-W) . (5.17) 
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This formula is the same given in [28] but it is different from the formula of Ahn 
and Nepomechie [32] in which L (9) is used instead of L{6). In Figure [2] log g a is 
plotted against log™. 

In the UV limit 0b — > — oo, the integrand is non-vanishing if 6 — > — oo and 
similarly, in the IR limit #b — > oo the non-zero contribution comes from the 6 — > oo 
domain. Using the fact that d#$(6>) = — 27r we find from the TBA equations 
(EH)) that 



L(-oo) = log(l + V2) , L(oo) = 0, 
Lq(-oo) = log(2 + V2) , L (oo)=log2. 



(5.18a) 
(5.18b) 



Now if one of the boundaries has reflection factor Rq and the other one has Rj, then 
from (IX25D . (OP and (EH)) (0) = -$(0 - B ), so 



\Pa / Vfl'1,1 



(5.19) 



since the factor gp = g^i cancels. Thus the TBA predicts the flow (1,3) — > (1,1), 
which is different from the prediction of [12]. We will see that TCSA supports our 
result. 



-0.2 



-0.4 



-0.6 



-0 . 8 



-6 



-2 



Figure 2: TBA: log(g) along the flow 



6 Flows with TCSA 

TCSA can be used to explore the fixed point of the renormalisation group flows. 
These flows can be implemented by varying the volume r while keeping the coupling 
constant K fixed in the Hamiltonian <\4M . Equivalently - and we choose this way 
- one can keep r fixed and vary k on some interval. Starting from k — 0, which is 
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the ultraviolet (UV) limit, the matrix h can be diagonalised at different values of k. 
The flow approaches a fixed point, that is a new supersymmetric conformal bound- 
ary condition. We should observe the eigenstates rearranging themselves into some 
degeneracy pattern from which one can identify the modules (the boundary condi- 
tions) using the characters and weight differences of the supersymmetric minimal 
model. 

It is important that the errors of the TCSA diagonalisation cannot be controlled 
easily. For example, it may happen that before the flow reaches the scaling region 
the truncation errors start to dominate. If we use various cuts and find that the 
flow picture does not change drastically (only the precision of the result gets higher 
with higher cuts), then it means that the unpleasant case mentioned above does not 
happen. 

Of course one can not establish the endpoint of the renormalisation group flow 
using TCSA. What one can see is that the flow goes in the vicinity of some super- 
conformal boundary condition. The exact infrared fixed point can never be reached 
by TCSA because of the truncation. We are looking for the range where the TCSA 
trajectory is closest to the fixed point. 

The characters of the (1,3) and (1,1) highest weight representations are 

X i, 3 {q) = 1 + q 1/2 + q + q 3/2 + 2q 2 + 2q 5/2 + 2q 3 + 3q 7/2 + V + 5g 9/2 + 5q 5 

+ 6g 11/2 + 8g 6 + 9g 13/2 + 10g 7 + . . . (6.1a) 

Xi,i(q) = 1 + q 3/2 + q 2 + <? 5/2 + q 3 + q" 2 + 2g 4 + 2q 9 ' 2 + 2g 5 

+ 3q 11/2 + 4g 6 + 4g 13/2 + Aq 7 + . . . (6.1b) 




(a) Integer levels (b) Half-integer levels 



Figure 3: The RG flows starting from the (1,3) NS and NS b. c.-s in SM(2,8) 

In the TCSA calculations the dimension of the Hilbert space was 393 for the 
integer levels and 344 for the half-integer ones. In Fig. [3] the normalised energy dif- 
ferences, (e~ e )/(ei— e ) are plotted, and the degeneracies are shown in parentheses. 
We have found that starting with the integer (half-integer) levels of the module (1,3) 
the flow tends to a degeneracy pattern corresponding to the integer (half-integer) 
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levels of the module (1,1). Since we did not perturb the (1,1) boundary we can 
assume that this boundary condition remained the same during the flow. Taking 
into account that (1,1) is bosonic and (1,3) fermionic, this - keeping in mind the 
conclusions after equation (12 . 91) - leads to the conclusion that the flows are 

(1,3) NS — »(1,1) S§J (6.2a) 
(1,3)^— >(1,1) NS (6.2b) 

in perfect agreement with the results of the TBA analysis based on our proposed 
reflection factor. 

The flow of the boundary entropy, log(g) can also be explored numerically using 
TCSA [TJ]. The trace in the partition function (12.91) can be approximated by the 
finite sum over the TCSA eigenvalues. For each value of / one must find the scaling 
region, the domain in r for which logZ behaves as in (15.161) . The volume r should 
be great enough for this scaling behaviour, but for too large values the truncation 
errors start to dominate and spoil this simple linear form. In Figure 4(a) it can 
be seen that the scaling region is centered around r/l = 3 and in Figure 4(b) the 
logarithm of the partition function is plotted. 

The partition function will depend only on the combination x = /iL 1 / 2- ^' 3 = 
hL 3 ^ = k/ 3 / 4 . After making the linear fit in r/l along the scaling region for different 
values of x one gets the product of the (^-functions as a function of x. However this 
function generally contains a linear term in L which corresponds to the free energy 
density coming from the boundaries. In order to get the correct final result this 
term should be subtracted. We extracted this term numerically from the large L 
behaviour of the naive fits for the ^-function. 

The numeric results are compared with the TBA data in Table [2] and they 
are plotted in Figure [5] together. Unfortunately there is no systematic method for 
determining the TCSA errors, but one can make estimates for them. One way is 
comparing the results obtained at different cuts, i. e. at different dimensions of the 
Hilbert space, here we found that for x < 1 the relative error is 0.5-2%. Another 
source of error is in the choice of the scaling region, the corresponding error is 0.2- 
0.7%. Finally, the error of the linear fit itself is about 0.2-0.4%. So our estimate 
for the relative error in the TCSA results is about 1-3% and the TBA and TCSA 
values for the ^-function agree within this error. Apart from the difference due to 
the truncation errors in gxcsA at large values of x (as in [H]) there is an excellent 
agreement between the two approaches. 
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-0.85 




(a) d(logZ)/d(r/Z) (b) logZ vs. r/l 



Figure 4: The partition function (x = 0.1) 



log a; 


logfi-TBA 


logfi-TCSA 


-4.605 


-0.0900 


-0.0932 


-3.912 


-0.1013 


-0.1059 


-3.507 


-0.1122 


-0.1203 


-2.996 


-0.1337 


-0.1449 


-2.302 


-0.1846 


-0.2008 


-1.609 


-0.2768 


-0.2941 


-0.916 


-0.4264 


-0.4322 


-0.693 


-0.4864 


-0.4906 


-0.224 


-0.6215 


-0.6177 





-0.6837 


-0.6922 


0.405 


-0.7817 


-0.7457 


0.693 


-0.8347 


-0.9015 


0.916 


-0.8659 


-0.9377 


1.504 


-0.9169 


-1.2646 



Table 2: TBA and TCSA results for logg 



7 Generalisation to SM(2, An + 4) 

The supersymmetric Lee- Yang model is the first member of the series of super- 
conformal minimal models with p = 2, the models SM(2,4n+4). Using TCSA we 
have found that for k > starting from the even (odd) levels of any NS module we 
end up at the even (odd) levels of the (1,1) module (see Figures El [7j). So the IR 
fixed point seems to be always the (1,1) boundary condition and depending on the 
fermion parity of the field corresponding to the UV boundary condition the flow is 
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-4 -3-2-10 1 2 



Figure 5: log(g) vs. log(x) along the flow, TCSA (dots) and TBA (solid line) results 
of type NS^NS (NS -> NS) or NS^ NS (NS ^NS): 



(1,3)ns- 


^)ns ' 


(7.1a) 


(1, 3)g§ - 


-(M)ns, 


(7.1b) 


(l,5) N s - 


->(M)ns, 


(7.2a) 


(!.5)ns- 


-(M)nsp 


(7.2b) 


(1,7) ns - 


(■*■' -0ns ' 


(7.3a) 




-(M)ns 


(7.3b) 




2 4 6 8 10 12 14 2 4 6 8 10 12 14 

(a) Even levels (b) Odd levels 

Figure 6: The RG flows starting from b. c. (1,5) in SM(2,12) 
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1 



1 



2468 10 2468 10 

(a) Even levels (b) Odd levels 

Figure 7: The RG flows starting from b. c. (1,7) in SM(2,16) 

8 Conclusions 

In this paper we studied the supersymmetric Lee -Yang model (SM(2,8) or 
M(3,8)) in the presence of boundaries, on a cylinder. First we proposed reflection 
factors for the NS type boundary conditions by considering our model as a reduction 
of the supersymmetric sine-Gordon model. The reflection factors of the SSG model 
contain two boundary parameters but a consistent reduction requires a functional 
relation between them, leaving only one boundary parameter which is related to the 
coupling constant of the boundary perturbation of the SLY model. In the massless 
limit even this boundary parameter gets eliminated. 

After determining the massless limit of the reflection factors we compared the 
energy levels predicted by the Bethe-Yang equation with the numerical spectrum 
calculated with the Truncated Conformal Space Approach. We have found very good 
agreement for the even levels of the super Verma module. However, we could not 
find any reflection factor that could describe the behaviour of the odd energy levels. 

Then we turned to the question of the boundary renormalisation group flows of 
the SLY model. The (1,3) boundary was perturbed by the field (5-1/201,3 which is a 
relevant, integrable perturbation that preserves supersymmetry. Using our reflection 
factor we wrote down the boundary TBA equations and determined the flow of the 
(^-function or equivalently, the boundary entropy. Comparing the UV and IR values 
of the (/-function we concluded that the boundary condition (1,3) flows to the (1,1) 
boundary condition. 

We checked this result using TCSA, diagonalising the Hamilton-operator of the 
system at various boundary couplings. From the state content and degeneracy pat- 
tern of the Hilbert space we could identify the final boundary condition as the 
boundary condition (1,1), in perfect agreement with the prediction of the TBA 
analysis. Finally we calculated the ^-function along the flow with TCSA, by approx- 
imating the partition function and extracting the (/-function from the scaling region. 
The TBA and TCSA results for the change of the boundary entropy along the flow 
are in excellent agreement. 
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The proposed reflection factors are different from those of Ahn and Nepomechie 
[T2] but the energy spectrum, the fixed point of the boundary renormalisation group 
flow and the TBA (/-function based on them are in very good agreement with the 
TCSA results. 

At the end of the paper we examined the boundary flows in the generalisations 
of the SLY model, the superconformal minimal models SM(2,4n + 4) with TCSA. 
We found that every Neveu-Schwarz (NS and NS) boundary condition flows to the 
(1,1) boundary condition. 
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A The SSG S-matrix 

The nth breather in the supersymmetric sine-Gordon model has mass 

/n7r\ 

m n = 2msmy— j. (A.l) 
The SSG breather S-matrix can be written in the form 

S^(e)=sM{6)SjM y {6). (A.2) 

Here [H] 

c(M')/m sinh(^) + isin(^7r)sinh(^) +isin(9f7r) 



sinh(#) — % sin(^j-7r) sinh(0) — i sin(^-7r 



^ sin^^ + ^cos^^ + zf) 

is the sine-Gordon S-matrix and the supersymmetric factor is 

S^ Y (6) = M^(6)G^\6) (A.4) 

with 

M^ m \e) = 

/ sin(^) + sin(^) \ 

1+1 sinh(e) U U coBh(f) 

sinh(e) ' sinh(|) 

o i vggM^lj 1+ .sin (S )-y^) 

smh(|) sinh(0) 

{ v^m^m o o -i+i ^m^n 

\ cosh(|) smh(0) / 

(A.5) 
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and 



f n+m \ l\ _ n— m ' 
G (n, m) ( m = 41 4A ^ 2 / A 6) 

9(2) 



sinh(|) I /• dt sinh(At) sinh((l - A)t) /itfl\ 

sinh(f) + zsin(Avr) 6XP I 7 t cosh 2 (§ ) cosh(t) Sm \ vr / 







(A.7) 



B SSG reflection factors 

The sine-Gordon reflection factor for the nth breather is Q29J) 

R^(e) = 4 l \e)R^\e), (B.i) 

where 

R (n) {9) = ( _ ir+ iCos(| + f ) cos(| - I - f ) sin(! + f) 



cos(|-f)cos(| + f + f)sin(|-f; 



^ sin(f + ^)cos 2 (|-f-^) 
llsin(f-||)cos 2 (| + f + S)- ( • j 



R ( ( l \9), which contains the boundary parameters r\ and $ is different depending on 
whether n is even or odd: 

Rf n \9) = S {2n \r 1 ,9)S {2n \id,9) , (B.3) 

where 

S ( 2n) / n\ _ TT sin(f)-cos(f-(/-|)f) sin(f)-cos(f + (/-|)f) 

1 ' J sin(f ) + cos(f - (/ - |)f ) sin(f ) + cos(f + (/ - |)f ) 1 " J 

and 

Rf n ~ 1 \e) = S^-^ir), 9)S {2n - l \ti9, 9) (B.5) 



with 



^n-i), ^ _ cos(f)-sin(f) !p3 s in(f)-cos(f-^) sin(f ) - cos(f + g 



cos(f ) + sin(f ) xi sin( |) + 006(2 _ M) sin (|) + C0S (| + 

(B.6) 

In the supersymmetric sine-Gordon model the reflection factors are of the form 

mm) 

R^ G (9) = R$l(9) ® R S vM9) . (B.7) 
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In the so-called BSSG + case, in which the fermion parity is conserved during the 
reflection 

'A + 



iW(0)=(Y rj (B.8) 
with 

A ± (9) =p(9 + i^)p(e- ^ P -) V2K(26)2-^ cos T f ) , (B.9) 
where for the fcth breather pk = n — kj. The functions K(9), P(6) are defined as 

A7g) = -1- TT F(A: ~ 1 + ^ )r(fc ~ ^ (BIO) 

00 r 2 ('A---M 

P (#) = TT^ 1 7i ; 1 r^/i 6 ^- }- (B.ll) 

k = l V 4 V 4 27TI / 

They have the following integral representations: 



00 



™-^!hJ$£bs>- (B ' 12) 







1 [dt smh(&) 
4 / t cosh 2 '*' 



= / - exp < - 7 / -— T27T7 > • (B.13) 

4> 



'cosh(f) k (| 
Then _Rsusy(#) can be written as 



yjgj 1 /-dt cosh(f) + cosh 2 (|) /ft 

Psusy{9) - ; ■ = exp <; - 7 / — — — jtj— — jtj— smh 



^/cosh(0) 1 4J t cosh 2 (f)cosh 



-2' 



/TT 



x r ( |~ 5> cos( | + }) ) ■ (B.14) 
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